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Chiral measurements with the Fixed-Point Dirac operator and 
construction of chiral currents 

P. Hasenfratz a * , S. Hauswirth a , K. Holland 3 -, T. J6rg a , F. Niedermayer a 

a Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, CH-3012 Bern, Switzerland 

In this preliminary study, we examine the chiral properties of the parametrized Fixed-Point Dirac operator 
D FP , see how to improve its chirality via the Overlap construction, measure the renormalized quark condensate 
£ and the topological susceptibility Xt, and investigate local chirality of near zero modes of the Dirac operator. 
We also give a general construction of chiral currents and densities for chiral lattice actions. 



1. Introduction 

There has been enormous recent interest in lat- 
tice QCD with chiral fermions, with several ap- 
proaches now being used in simulations |Q. We 
use one approach, Fixed- Point actions ||, which 
have many desirable features, including chiral 
symmetry. This preliminary study gives the first 
chiral measurements with this action. We also 
show how to construct chiral densities and con- 
served currents for chiral actions. 

2. Chiral properties of D FF 

A lattice Dirac operator D which satisfies the 
Ginsparg- Wilson (GW) relation || 

{D-\ l5 } = 2aR l5 , (1) 

has exact chiral symmetry at non-zero lattice 
spacing a Bj. The exact Fixed-Point Dirac op- 
erator satisfies exactly the GW relation || with 
a non-trivial local function R. The function 
2R can be absorbed into the definition of D 
(6) and, for notational simplicity, we shall write 
2R = 1 in most of the following equations. The 
parametrized Fixed-Point Dirac operator D FP , as 
described in Jjj , is an approximate solution of the 
GW relation. In Fig. |], we show the eigenvalues 



of D for Fixed-Point gauge action S g 
figurations || of volume 4 4 and 5 4 with lattice 
spacing a « 0.16 fm (as measured by the Som- 
mer parameter tq). For an exact solution D of 
the GW relation, the complex eigenvalues of D 

*Based on a talk by T.J. and a poster by K.H. 



lie on the circle of radius 1 and center at (1,0). 
As the eigenvalues of D FF lie very close to the 
GW circle, we see that D FP satisfies the GW re- 
lation well. 

The GW relation is equivalent to 

A*A=l, A = l-D. (2) 

In Fig. |[ we plot the 10 smallest eigenvalues of 
A^A divided by the largest eigenvalue of A' A for 
10 Fixed-Point gauge configurations with lattice 
spacing a k 0.13 fm. The largest eigenvalue is ~ 
1.5 and 38 for L> FP and D wilson respectively. We 
see that the A' A eigenvalues are much closer to 1 
using D FP than by using £) Wllson . For larger vol- 
umes, A^ A is more likely to have a few very small 
eigenvalues, but the vast bulk of the eigenvalues 
is close to 1. 

An exact solution of the GW relation is given 
by the Overlap operator || 

^Overlap = 1 _ A /^JTa, A=1-Dq. (3) 

The already very good chiral behavior of D FF can 
be made exact by using the Overlap construc- 
tion. Very precise chiral symmetry is required for 
some measurements, such as the chiral conden- 
sate. Using D FF , A* A is close to 1 and we use 
a polynomial approximation of 1/V A^ A, which 
should converge rapidly. In Fig. ||, we see that 
the breaking of the GW relation, given by 

B 2 = \(D + D ] - D^D).v\ 2 , 



(4) 



where v is a random vector, falls off exponentially 
as we increase the order of the polynomial approx- 
imation of 1/VaTA and becomes < O(10 -10 ) at 
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Figure 1. Eigenvalue spectrum of D . 



Figure 2. Chirality breaking of D 
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polynomial order 10 for this volume. If the eigen- 
values A of Z)Ovoriap lie exactly on the GW circle, 
then A = A/(l — A/2) is purely imaginary. In 
Fig. |], we see that real(A) falls off exponentially 
as we increase the polynomial order, as the eigen- 
values A get closer and closer to the GW circle. 
The most rapid decrease is at the finest lattice 
spacing a « 0.10 fm. 

A lattice action must be local for the continuum 
results to be universal i.e. independent of the de- 
tails of the lattice action. Locality means that 
fields at large separation have an exponentially 
small coupling. Optimizing the locality is essen- 
tial so that e.g. the exponential fall-off of correla- 
tion functions can be separated from direct cou- 
plings in the lattice action. The locality of a lat- 
tice Dirac operator can be measured by 



f(r) =max {\D.v\, \\y - x\\ = r}, 



(5) 



where v is a vector with point source at x and 
r is the square norm. In Fig. [5], we see that 
the Overlap operator constructed from the inter- 
acting D FP is more local than by using the free 
Z) Wllson , which in turn is more local than by using 
the interacting D Wilson . 



3. Chiral condensate 

For a theory of Nf > 2 massless quark flavors, 
the chiral symmetry is spontaneously broken if 
the chiral condensate {ipip) ^ 0. Using chiral per- 
turbation theory or random matrix theory, the 
chiral condensate at finite quark mass and vol- 
ume has been calculated in the continuum, both 
for quenched and full QCD jlo|]. The quenched 
QCD condensate E mj y,Q = — (^ip}m,v,Q in topo- 
logical sector Q at small m is 



E m ,v,Q = \Q\/mV + 

mVY?[I M (z)K M {z) + I lQ \ +1 (z)K ]Qhl (z)], 



(6) 



where Iq and Kq are modified Bessel functions, 
z = rnEV and E = —{ipt/j) is the quantity we 
want to measure. The term \Q\/mV comes from 
the zero modes, whose contribution is not sup- 
pressed by the fermion determinant in quenched 
QCD at finite volume. Subtracting the contribu- 
tion of the zero modes, for Q ^ 



lim 



y sub 



(7) 



By measuring this ratio for several masses, vol- 
umes and topological sectors, the scaling behavior 
can be used to extract E lllfl . 
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Figure 3. Breaking of Ginsparg- Wilson relation. Figure 4. Deviation of £)° vcrla P eigenvalues from 

GW circle. 



The subtracted condensate is measured by 

V£y,Q = Tr sub [( J D(m)2 J R)- 1 - 1/2] 

D(m)2R = (1 - m/2)D(0)2R + m, (8) 

where D and R are related via the GW rela- 
tion. We use the Overlap construction with D FP 
(hence the Overlap D FP ) so that the chiral sym- 
metry is precise enough to go to quark mass 
ma = 1CP 4 . The trace is measured stochasti- 
cally using random Z(2) vectors. The topology 
Q is determined from the chirality of the lowest 
eigenmodes of D^D, found with an Arnoldi solver. 
For the topology scan, a Legendre polynomial of 
order 2 is used to approximate X/\/ A>A. As a 
byproduct, we determine the topological suscep- 
tibility from an ensemble of 200 10 4 configura- 
tions to be r\\t — 0.0612(75) corresponding to 
Xt = (196 ± 6 MeV) 4 . In the volumes we exam- 
ine, all Q zero modes have the same chirality and 
their contribution to the condensate is subtracted 
by using random Z(2) vectors whose chirality is 
opposite to those of the zero modes. 

We measure the condensate in volumes 6 4 ,8 4 
and 10 4 in topological sectors \Q\ = 1,2 at lat- 
tice spacing aw 0.13 fm. The statistics are given 
in Table 1. We use 10 random Z(2) vectors to 
measure the trace for each configuration and a 
multi-mass solver to invert D simultaneously at a 



number of masses. We use Legendre polynomials 
of order 5, 7 and 10 to approximate X/V At A for 
volumes 6 4 ,8 4 and 10 respectively. This gives 
sufficiently precise chiral symmetry — increasing 
the polynomial order further, the relative change 
in S™ b yQ /m is < O(10" 4 ). The 10 smallest A 1 A 
eigenvalues are projected out and treated exactly. 
In Fig. ||, we plot a 3 E™ b l/ g/ma as a function 
of ma for several volumes and topological sec- 
tors. We see that the data reach a plateau at 
small quark masses. Including lattice artifacts, 
the finite-size scaling behavior is 

^™y,Q ci c 2 m 

m 2\Q\ a 2 a ' ( ' 

where ci,C2,... are unknown coefficients which 
have to be fitted. It is natural to assume that 
these coefficients, dependent on the UV fluctua- 
tions, are independent of the topological charge 
Q. As shown in Fig. 0, we fit the data at ma = 
10 -4 and determine that the bare condensate is 
a 3 T, = 4.68(26) x 10~ 3 . We can convert this us- 
ing the Sommer parameter (ra/a = 3.943(60) has 
previously been measured at this bare coupling), 
which gives r^E = 0.287(16)(13), where the first 
error is statistical and the second the uncertainty 
in the scale. 
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Figure 5. Locality of D° vclla v using £> wilson and 
D FF (exponential fall-off with v ~ 1.6). 



The renormalization factor Z5 for the chiral 
condensate is related to the quark mass renormal- 
ization by Z$ = l/Z m . The mass renormaliza- 
tion can be determined from hadron spectroscopy 
measurements, as done in fl^| . Hadron mass mea- 
surements using the Overlap D FF have been done 
at the same lattice spacing a w 0.13 fm giv- 
ing Zs — 0.80(6) in the renormalization group in- 
variant scheme. Using this, the renormalized chi- 
ral condensate is rgE = 0.230(13)(10)(17), with 
statistical, scale and renormalization errors re- 
spectively. This gives £ = (242 ± 9 MeV) 3 , or in 
the MS scheme, £^(2 GeV) = (270 ± 10 MeV) 3 . 
Comparing our measurement with other recent 
determinations of the renormalized chiral conden- 
sate [ p^JT^ I, we see in Fig. || that there is good 
agreement. 



4. Local chirality of near zero modes 

From instanton physics, it is expected that 
near zero modes of a lattice Dirac operator are 
localized around instantons and anti-instantons 
and that, in these regions, the modes are pre- 
dominantly cither left- or right-handed. These 
near zero modes form a finite density p(0) as 
V — ► 00, breaking the chiral symmetry via the 
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Table 1 

Statistics for the condensate measurement. 
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Figure 6. a Y, m y \Q\/ma vs. ma for several vol- 
umes and topological sectors. 



Banks-Casher relation (ipip) = — np(0). A mea- 
sure of the local chirality at lattice site x is |H| 



tan -(1 



X(x))} 



IpR'tpB.(x) 



(10) 



If the near zero modes are localized and locally 
chiral, then where ip'ip(x) is large, X should be 
close to ±1. If the modes are not locally chiral, 
then X should be close to 0. 

We analyzed the 10 smallest near zero modes of 
the Overlap D FF for 60 10 4 configurations at lat- 
tice spacing a sw 0.13 fm, using a Leg endre poly- 
nomial of order 2 to approximate 1/v A^A (these 
modes were already used to determine the topol- 
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Figure 7. a?Y, mV _\Q\/ma at ma = 10 4 vs. 
V/(2\Q\a i ) and fit to finite -size scaling behavior. 



Figure 8. Comparison of measurements of 
roSj^g(2 GeV) by different groups. 



ogy Q of the gauge configurations). In Fig. |], we 
plot the distribution P(X), where we include the 
top 1%,5% and 10% lattice sites x where i/j'i/j(x) 
is largest. We see a clear double-peaked distri- 
bution, with peaks far from X — 0. Where the 
modes are most localized, they are also very chi- 
ral. This supports the picture of instanton domi- 
nance of the near zero modes, in agreement with 
16 1, unlike the original claim of |15|. 



5. Construction of chiral currents 

For a chiral lattice action, the advantages of us- 
ing the conserved currents are that Zy = Za = 1, 
the current and densities arc automatically 0(a) 
improved and expressions have the same form as 
in the continuum. The axial-vector current with 
the Overlap operator was already given in Jl7| . 
Here we define conserved currents (obtained via 
the Noether procedure) which are easy to use 
in numerical simulations. The vector and axial- 
vector currents are 

V£(x) = ^r a [P L K^x)P R + P R K^x)P L }^ 
A«(x) = ijT a [P L K^x)P R -P R K^x)P L ]^ 
P RtL = (1 ± 75 )/2 P L , H = P L , R ± 75 D/2. (11) 



These currents are conserved and have the correct 
properties under global vector and axial-vector 
transformations e.g. 



*-?—Vhx) = if abc Al{x) axial-vector 



(12) 



The kernel K^iXjU) which appears in the cur- 
rents is defined by 



Kp(x,U) 



.5D{U»{x)) 
Sa fi (x) 



\a=0 



(13) 



The kernel can be calculated easily — the gauge 
link Ua(x) receives an infinitesimal U(l) rotation 
and the Dirac operator is calculated on this new 
and original gauge configuration (differing only 
at one link). The ratio of the difference of these 
two Dirac operators to the U(l) variation gives 
K[i(x,U). Details of the derivation and proper- 
ties of the currents and densities will be given in 
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Figure 9. Distribution P(X) for top 1%,5%,10% 
lattice sites x with largest ip'tp(x). 



6. Summary 

The parametrized Fixed-Point Dirac operator 
D FP has very good chiral properties and low- 
order polynomial approximations can be used in 
the Overlap construction to make the chiral sym- 
metry very precise. The Overlap operator is 
also much more local using D FP than by using 
£)Wiison_ Tj sm g finite-size scaling, the bare chiral 
condensate is measured as rg£ = 0.287(16) (13) 
and, estimating the renormalization factor from 
hadron spectroscopy, the renormalization group 
invariant condensate is r^S = 0.230(13)(10)(17), 
i.e. S = (242±9 MeV) 3 . In addition, the topolog- 
ical susceptibility is Xt — (196 ± 6 MeV) 4 . Near 
zero modes of the Dirac operator are locally chi- 
ral, which supports the picture that they are lo- 
calized on instantons. We give a general and 
practical construction of chiral currents and den- 
sities, to be used in e.g. determining renormaliza- 
tion factors and decay constants. 
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